Robust ultrafast currents in molecular wires through Stark shifts 



O 

o 



Ignacio Franco,^ Moshe Shapiro,^ and Paul Brumer^ 

^Chemical Physics Theory Group, Department of Chemistry, 
and Center for Quantum Information and Quantum Control, 
University of Toronto, Toronto, Ontario, Canada. 
^Chemical Physics Department, The Weizmann Institute, 
Rehovot, Israel, and Depts. of Chemistry and Physics, 
The University of British Columbia, Vancouver, B.C., Canada 
(Dated: February 1, 2008) 

A novel way to induce ultrafast currents in molecular wires using two incident laser frequencies, 
u) and 2uj, is demonstrated. The mechanism relies on Stark shifts, instead of photon absorption, to 
transfer population to the excited states and exploits the temporal profile of the field to generate 
phase controllable transport. Calculations in a frans-polyacetylene oligomer coupled to metallic 
leads indicate that the mechanism is highly efficient and robust to ultrafast electronic dephasing 
processes induced by vibronic couplings. 
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Considerable effort has been devoted to studies of the 
properties of molecular wires [H, as the basic building 
blocks of nanoscale electronics. The focus is normally 
placed on the transport properties of metal-molecule- 
metal junctions subject to a bias voltage. In this regime, 
the metallic leads are the main source of electrons for the 
transport while the molecular system serves as a trans- 
porting medium that can be chemically functionalized to 
modify the I-V characteristics of the junction. 

Here we consider an alternative situation in which 
the junction is not subject to a bias voltage and where 
the molecule serves as the main source of transporting 
electrons. The composite system is taken to be spa- 
tially symmetric and net currents are induced using lasers 
with frequency components w and 2lu. Such fields give 
rise to phase-controllable transport in symmetric sys- 
tems 0, S 3 even when they have a zero temporal mean. 
This rectification effect first appears in the third order re- 
sponse of the system to the incident radiation The 
setup is of interest in molecular wires since it can pro- 
duce ultrafast currents and could lead to the develop- 
ment of molecular switches that operate on a femtosec- 
ond timescale. 

In the usual way that the scenario is applied the 
frequencies of the uj + 2uj field are chosen at or near 
resonance. The resulting multiphoton absorption pro- 
cesses serve the double purpose of creating laser-induced 
rectification and transferring population to transport- 
ing states. Such setup, however, is fragile to decoher- 
ence processes since it relies on creating coherent su- 
perposition states Q. This aspect becomes particularly 
troublesome when applying the scenario to molecular 
nanojunctions as most molecules exhibit strong elcctron- 
vibrational couplings that introduce ultrafast coher- 
ence loss (in less than 10 fs 0]) and internal relaxation 
mechanisms in the electronic dynamics 

In fact, recent calculations f7l| show that vibronic cou- 



plings can make the laser rectification in molecular wires 
extremely inefficient. In that study we found that less 
than 4% of the photoexcited electrons participated in the 
current, whereas equivalent but rigid systems exhibited 
60% efficiencies. Hence, in order to generate sizable cur- 
rents considerable energy from the laser field needs to be 
dumped into the nanojunction compromising its struc- 
tural integrity. Further, employing faster or stronger 
pulses is not very helpful Q since they introduce other 
satellite channels in the control and do not appreciably 
overcome the deleterious effects of the vibronic couplings. 

In this Letter we present an alternative mechanism 
that is remarkably robust to electron-vibrational cou- 
plings, survives in the presence of decoherence and ther- 
mal effects, and is able to induce large currents in molec- 
ular wires with efficiencies > 90%. Instead of promoting 
electrons to the conduction band through absorption of 
photons, we work far off resonance and employ Stark 
shifts to nonadiabatically couple the ground and excited 
electronic states. Phase-controllable symmetry breaking 
is achieved by exploiting the difference in the intensity of 
io + 2lu fields for positive and negative amplitudes. 

As a system we consider a trans-polyacetylene (PA) 
oligomer connected by its ends to macroscopic metallic 
leads. The leads are treated as rigid semi-infinite tight 
binding chains while the oligomer is described using the 
well-known Su-Schrieffer-Heeger (SSH) Hamiltonian 
The SSH model treats the molecule as a tight-binding 
chain in which the tt electrons are coupled to distortions 
in the polymer backbone by electron-vibrational inter- 
actions. It is remarkably successful in reproducing the 
electronic structure and dynamics of excitations of PA. 

The photoinduced electron-vibrational dynamics of 
the nanojunction is followed in a mean-field (Ehren- 
fest) mixed quantum-classical approximation where tran- 
sitions between instantaneous eigenstates are allowed [§] . 
The effect of lattice fluctuations is incorporated by propa- 
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gating an ensemble of quantum-classical trajectories with 
initial conditions selected from the ground-state nuclear 
Wigncr distribution function of the chain. In this way 
the dynamics reflects the initial nuclear quantum distri- 
bution and is subject to the level broadening and internal 
relaxation mechanism introduced by the vibronic cou- 
plings in the wire |71]. We have found that this method 
describes the trends observed in a recent sub-10 fs exper- 
iment performed on PA [lo| . In addition, we introduce 
projective lead-molecule couplings so that basic aspects 
of the Fermi sea are taken into account without explicitly 
following the dynamics of the essentially infinite number 
of particles in the leads. 

The nanojunction is described as a one-dimensional 
lattice in which each site n corresponds to the position 
of an atom, and is defined by the Hamiltonian 



Hit)^Hi^+Hs.L + Hs{t)+Hs.R + HR. (1) 



Here Hi^ 
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H.c. and Hn 



H.C. describe the left (L) and 
right (R) lead with hopping parameter ^load- The oper- 
ator cjj J, {cn,s) creates (annihilates) a fermion in site n 
and spin s, and H.c. denotes the Hermitian conjugate. 
The SSH oligomer with Hamiltonian Hs{t) — Hc\ -l-i?iatt 
is situated between sites n = 1, • • • , iV and is coupled to 
an electric field E{t) in the dipole approximation. Specif- 
ically, the electronic part of -ffs is 



JV-l 



n— l,s 

N 



(2) 



where to is the hopping integral for zero displacement, w„ 
is the monomer displacement of site n and a describes 
the electron-ion coupling between neighboring sites. In 
turn, x„ = (na + Un) is the position operator for site n, a 
the lattice constant, and — |e| the unit electronic charge. 
The lattice is described by 
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with force constant K and (CH) group mass M. 

We consider the case where no bias voltage is applied 
across the bridge and focus on the short time dynamics 
of the system. In this regime, the basic function of the 
leads is to absorb electrons from the molecule with en- 
ergy higher than the Fermi energy ep, taken to be the 
zero-reference energy, and to block them otherwise. We 
focus on electronic transport and model this effective cou- 
pling using projection operators. Specifically, the lead- 



molecule coupling 

Hs-L{t) = -icoup ^?,«(04,.c„.. +H. 



(4) 
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is a restricted tight-binding interaction with coupling 
constant icoup in which only electrons with energy 
e > ep are deposited in the contacts. Here V^{t) = 
>eF lTs(i))(7s(OI projects into the instantaneous 
molecular tt* light-dressed eigenorbitals defined by the 
eigenvalue problem i?oi(t)|7s(i)) = e-f\ls{t)), andP^ ^ = 
{n\P^\m) with \n) = cJj^jO) where |0) is the vacuum 
state. Note that since the SSH model is electron-hole 
symmetric a subsequent inclusion of hole transport is just 
expected to double the resulting current. 

In the mean-field approximation Q the nuclei are de- 
scribed by classical trajectories determined by 



+ 2aRe {pn,„+i - p,u„-i} - \e\Eit) (p„^„ - 1) . (5) 

Here p„^„ = Ee.s is the elec- 
tronic reduced density matrix, and /(e, s) is the time- 
independent initial distribution function (that takes val- 
ues or 1 depending on the occupation of each level with 
energy e and spin s). We take the chain to be clamped, 
so that ui{t) = UN{t) = at all times. In turn, the or- 
bitals |e(t)) that form the A/'-electron wavefunction |5'(t)) 
satisfy the time-dependent Schrodinger equation 
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where |ep) = X^nsP^'^l^)!'^) describes the part of the 
orbital in region P=L, R or S. We take the leads and 
molecule to be initially detached and, since the projective 
coupling [Eq. (jl])] already incorporates the Fermi block- 
ade imposed by the leads, only consider orbitals that are 
initially in the molecular region, so that |e^(0)) = with 
(3 =L or R. 

As the composite system is spatially unbounded, it is 
necessary to obtain an equation of motion for the system 
part of the orbitals jeg) where we do not need to explicitly 
follow its behavior in the lead regions. This can be ac- 
complished by applying Laplace transform techniques to 
the evolution equation for the lead components, to ob- 
taim |e^(t)) = -i-/Jc//3(t _r)iJs-/?(r)|es(r))dr. Here 
the evolution operator for the isolated leads U^{t) — 
exp{—iHist/H) acts as the memory kernel of the con- 
volution integral. Its matrix elements in site repre- 
sentation can be determined in closed form. For this, 
note that H/^ is diagonalized by the basis transforma- 
tion Cn,s = (xt)^^^ Y^kLi sin [(" - nf3)kn/M] Ck,s, where 
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A4 is the number of lead sites, riL = 1 and tt-r 



TV. 



It then follows that U!^^{t) 



^ 1 



«"+™-'"^ J„+m-2„, (t^O for large M, where J„(z) = 
^IT" So 6^^'^°^^ cos{n9) dd is a Bessel function of the first 
kind of order n. Hence, 
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V / /C(i-r)r„,™(t,T)(m|es(T))dr, (6) 
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where r„,„(i,T) = Vl^{t)Vl^(T) + Vl^^{t)V%,jT). 
Equation ^ is an effective non-Markovian Schrodinger 
equation for the system part of the orbitals. The first 
term corresponds to the single-particle Schrodinger equa- 
tion for the isolated system. The convolution integral 
with memory kernel lC{t) = UQ Q{t) — t^jv+i jv+i(^) = 
2Ji (2ttoid^j/2t^^ describes the transfer of population 
from the system into the leads. 

We now invoke the wide bandwidth approximation, 
where ?i/2iioad becomes the fastest time scale in the prob- 
lem. In this regime, Eq. ([6]) reduces to its Markovian 
limit 
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and the leads are effectively mapped into a negative imag- 
inary (absorbing) potential. The projection operators in 
the coupling, contained within r^^(t) = r„ ,„(t,t), en- 
sure that only those electrons with sufficient energy get 
absorbed with proper conservation of the antisymmetry 
principle. The equation is valid for times t ^ fi/^iead, 
but it can be used for all times by slowly turning on the 
lead- molecule interaction; a strategy that we adopt. Note 
that the field influences the dynamics directly by induc- 
ing transitions among states, and indirectly by modifying 
the wire-lead couphngs. Equations ([5]) and ([7]) constitute 
a closed set of N{N + 2) first-order differential equations 
that are integrated using a Runge-Kutta method of or- 
der eight. The projective term in the dynamics ^{t) 
is updated at every time step by diagonalizing Ha\(t). 
Throughout we use the standard SSH parameters [8|: 
a = 4.1 eV/A, K = 21 eV/A^, to = 2.5 eV, M = 1349.14 
eV fs^/A^ and a = 1.22 A. In turn, we take the molecule 
and leads to be weakly coupled with icoup/^iead = 0.1 eV. 

Electronic dephasing due to vibronic coupling is in- 
corporated by integrating the equations for an ensem- 
ble of 1000 initial conditions. The initial conditions are 
generated using the following strategy 0]: the starting 
optimal geometry is obtained by minimizing the total 
ground-state energy of the chain. Then, a normal mode 
analysis around this geometry is performed, yielding the 
ground state nuclear Wigner phase space distribution 
function pyv{u,p) in the harmonic approximation, where 




FIG. 1: Time dependence of (a) the electric field; (b) the 
current entering the left and right leads and; (c) the instanta- 
neous field-dressed orbital energies for states near the energy 
gap for a 100-site fiexible wire under the infiuence of the field 
in Eq. ((8]l with (p^uj — '24>lj = 0. Note the bursts of charge 
deposited in the right lead when the field bridges the energy 
gap. Here j'l is so small that is barely visible. 



u = {ui, • • • , ujv) and p = M(mi, • ■ • , un)- By impor- 
tance sampling this distribution an ensemble of lattice 
initial conditions {m*(0),p*(0)} is generated. The associ- 
ated initial values for the orbitals {|e*(0))} are obtained 
by diagonalizing iJci for each initial geometry of the lat- 
tice {m*(0)}. Each member i of the ensemble defines 
a quantum-classical trajectory (m*(0),p*(0), |\E'*(0))) 
{u^ (t) , (t) , and the set is used to obtain ensem- 

ble averages. We note that the resulting initial state is 
basically stationary under field-free evolution. 

The current entering into lead /3 is defined by jp = 
— \e\dpi3/dt, where = ^n^pPn,n is the number of 
electrons in lead (3. Using the same set of techniques 
employed to arrive at Eq. ([7]), one can obtain an equa- 
tion for j/3(t) that solely depends on molecular properties: 



2lo pulses manifests 



Any rectification generated by to 
as jL - 7^ 0. 

We study neutral 100-site chains initially in the ground 
electron-vibrational state. The geometry of the chain 
consists of a perfect alternation of double and single 
bonds. The electronic structure is composed of 50 dou- 
bly occupied valence tt orbitals and 50 empty tt* states, 
separated by an energy gap of 2 A = 1.3 eV. We follow 
the electron-vibrational dynamics of the wire under the 
infiuence of an -f- 2lo field of the form 



m 



, C0s(tjt -I- 0^) -f £2^^ C0s(2tjt + 02tj)- (8) 



The field is smoothly turned on and off in 100 fs and has 
constant amplitude for 400 fs. The frequency hio = 0.2A 
is chosen far off-resonance from the system's interband 
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FIG. 2: Phase dependence of (a) the net rectification — ^r, 
where g,? = jj^ jp it)dt and; (b) efficiency ?? = (gL - ga)/ (gL + 
gn) of the process when field ((8l is appUed to fiexible (black 
dots) and rigid (gray dots) f 00-site wires. Note the assistance 
of phonons in the rectification. 

transition frequencies so that Stark shifts, and not pho- 
ton absorption, dominate the dynamics. The field am- 
phtude used is e2uj = 6.1 x 10"'^ V with = 2e2uj, 
which corresponds to an intensity l2uj ~ 5 x 10® W cm~^. 

Figure [T] shows the field, the currents and electronic 
structure of the chain, averaged over all trajectories, 
when the relative phase of the pulse is 4>2uj — '^'Pui = 0. 
The single-particle spectrum displays considerable Stark 
shifts and the field effectively reduces the energy gap 
of the oligomer, causing frequent crossings between the 
ground and excited electronic states in individual trajec- 
tories. Since the wire's ground state is non-degenerate 
and of definite parity, the lowest order contributions to 
the Stark effect is quadratic in the field. Hence, when 
|i?(t)| is maximum the energy gap acquires its minimum 
value. At the crossing times population is transfered from 
the valence to the conduction band and bursts of charge 
are deposited in the leads. 

Note that for (l)2uj—2(j)^ — almost all excited electrons 
are deposited in the right contact only. Symmetry break- 
ing arises due to the difference in the maximum for 
positive and negative field amplitudes. Even when E(t) 
has a zero temporal mean, it consists of narrow peaks 
with large \E{t) \ for positive amplitudes, and shallow and 
broad features when E{t) is negative (Fig. [T^). For this 
reason, the Stark effect is only sufficiently strong to close 
the energy gap when the field has a positive amplitude. 
Thus, transfer of population to the conduction band and 
absorption of electrons by the leads always occurs when 
the laser is pointing at a particular and the same di- 
rection, in this way inducing directed transport in the 
system. 

The phenomenon depends intimately on the relative 
phase. For instance, for the case of (j)2ui — '24>ui = tt/2 (not 
shown) the electric field exhibits equal intensity for pos- 
itive and negative amplitudes. The energy gap is closed 
every half period, and since the field changes sign be- 
tween consecutive interband couplings, no net current is 
produced. 

Figure [5] shows the net difference in charge deposited in 



the left and right leads after the pulse is over for different 
laser phases, as well as the efficiency of the process. For 
comparison purposes the plot also includes results ob- 
tained in an equivalent but rigid system. The system is 
made rigid by increasing the mass of the (CH) groups by 
10^. We first note that the mechanism is robust to deco- 
herence effects due to coupling to the vibrational degrees 
of freedom as well as satellite contributions due to para- 
site multiphoton absorption. In fact, 90% of the excited 
electrons can participate of the net current. Further, the 
sign and magnitude of the effect can be manipulated by 
varying the laser phases. By changing the relative phase 
by TT the magnitude of the effect stays the same but the 
direction of the rectification is reversed. 

In the flexible wire the rectification exhibits an almost 
sinusoidal dependence on <^2w — 20^. By contrast, in 
rigid wires for certain range of phases no currents are 
induced since the maximum field amplitude is not large 
enough to couple valence and conduction bands. Hence, 
in this range the currents observed in the flexible wire 
are phonon-assisted, i.e. the level broadening introduced 
by the vibrations permit the nonadiabatic coupling. The 
currents generated in the flexible case are always larger 
than in the rigid example. However, in rigid wires the 
mechanism can exhibit perfect efficiencies. 

In conclusion, we have demonstrated a way to effi- 
ciently induce phase-controllable transport along molec- 
ular wires on a femtosecond timescale using uj + 2u; fields. 
The effect employs Stark shifts to transfer population to 
transporting states, and exploits the difference in the in- 
stantaneous field intensity for positive and negative am- 
plitudes to generate rectification. It is robust to electron- 
vibrational couplings and other decoherence mechanisms, 
and is best suited for long oligomers. 
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